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Abstract 

The Q 2 evolution of fragmentation function in non-equilibrium QCD may be necessary to study 
hadron formation from quark-gluon plasma at RHIC and LHC. In this paper we study DGLAP 
evolution equation of fragmentation function in non-equilibrium QCD by using Schwinger-Keldysh 
closed-time path integral formalism. For quarks and gluons with arbitrary non-equilibrium dis- 
tribution functions f q (p) and f g (p), we derive expressions for quark and gluon splitting functions 
in non-equilibrium QCD at leading order in coupling constant a s . We find that the quark and 
gluon splitting functions in QCD medium depend on non-equilibrium distribution functions f q (p) 
and f g {p)- We make a comparison of these splitting functions with that derived by Altarelli and 
Parisi in vacuum and find that the non-equilibrium QCD splitting function carries an extra factor 
[1 + fA(PA)] 2 x [1 + /b(pb)] 2 x [1 + fc{pc)] 2 where parton A emits a parton B. 
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I. INTRODUCTION 



RHIC and LHC heavy-ion colliders are the best facilities to study quark-gluon plasma 
in the laboratory. Since two nuclei travel almost at speed of light, the QCD matter formed 
at RHIC and LHC may be in non-equilibrium. In order to make meaningful comparison 
of the theory with the experimental data on hadron production, it may be necessary to 
study nonequilibrium-nonperturbative QCD at RHIC and LHC. This, however, is a difficult 
problem. 

Non-equilibrium quantum field theory can be studied by using Schwinger-Keldysh closed- 
time path (CTP) formalism However, implementing CTP in non-equilibrium at 
RHIC and LHC is a very difficult problem, especially due to the presence of gluons in 



one-loop resumed gluon propagator in 
4|- 



non-equilibrium and hadronization etc. Recent 
non-equilibrium in covariant gauge is derived in 

High pt hadron production at high energy e + e~ , ep and pp colliders is studied by using 
Collins-Soper fragmentation function . For a high px parton fragmenting to hadron, 
Collins-Soper derived an expression for the fragmentation function based on field theory and 
factorization properties in QCD at high energy \t\ . This fragmentation function is universal 
in the sense that, once its value is determined from one experiment it explains the data at 
other experiments. 

Recently we have derived parton to hadron fragmentation function in non-equilibrium 
QCD by using Schwinger-Keldysh closed-time path integral formalism [a]. This can be 
relevant at RHIC and LHC heavy-ion colliders to study hadron production from quark- 
gluon plasma. We have considered a high p? parton in medium at initial time r with 
arbitrary non-equilibrium (non- isotropic) distribution function f(p) fragmenting to hadron. 
The special case f(p) = -tw — corresponds to the finite temperature QCD in equilibrium. 

e T ±1 

We have found the following definition of the parton to hadron fragmentation function in 
non-equilibrium QCD by using closed-time path integral formalism [8(. For a quark (q) with 
arbitrary non-equilibrium distribution function f q (k) at initial time, the quark to hadron 
fragmentation function is given by 



D H/q (z,P T ) = — — ~ = *~r fdx '-^— ! t 1 t' ' 1 ' " /: 

1 1 



2z[l + f g (k))J ' (2vr 



W-l 



^tr Dirac ^tr color [ 7 + < in\il)(x ,x T ) $[x ,x T ] (J H (P + } T )a H (P + ,0 T ) $[0] ip(0)\in >\1) 
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where z (=tf) is the longitudinal momentum fraction of the hadron with respect to the 
parton and Pt is the transverse momentum of the hadron. For a gluon (g) with arbitrary 
non-equilibrium distribution function f g (k) at initial time, the gluon to hadron fragmentation 
function is given by 

I , Ad-2 r 



Dh/ 9 (z,Pt) 



2zk+ [1 + f g {k)\ 



(27T 



id— 1 



1 



-£[< in\F^{x-,x T ) $[x-,x T ] 4(P + ,0 r )a„(P + ,0 T ) $[0] F, 



fia 



in>\. 



(2) 



0=1 



In the above equations \in > is the initial state of the non-equilibrium quark (gluon) medium. 
The path ordered exponential 

f-0 



= V exp[ig f dXn- A a (x^ + n^X) T a ] 

J — oo 

is the Wilson line Q, |q] . 



(3) 



Eqs. (CQ) and (j2J) can be compared with the following definition of Collins-Soper fragmen- 
tation function in vacuum [5]: 



D H/q (z,P T ) 



1 

2z 



d x - - XT c ik+x-+iP T -x T /z 



(2vr; 



d-l 



^tr Dirac itr color [ 7 + < 0\ij(x-,x T ) $[x",x T ] 4(P + , T )a H (P+, T ) $[0] $(0)|0 >] (4) 



and 



Z,P T 



2zk+ 



" Xt c ik+x-+iP T -x T /z 



(2n) d - 1 



£[< 0\F^(x-,x T ) $[x-,x T ] a^ H (P+,0 T )a H (P + ,0 T ) $[0] F+(0)|0 >]. (5) 
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a=l 



Since the fragmentation function is a non-perturbative quantity, we do not have theoret- 
ical tools in QCD to calculate it yet. The normal procedure at high energy pp, ep and e + e~ 



colliders is to extract it at some initial momentum sea 



11 



scale n by using the DGLAP evolution equation flO . 

v-Q^ D i^w( z ) = J2 J z — p ij(-^) D 3^/4y)- 



e up and then evolve it to another 
3: 



(6) 



In the above equation Pij(z) is the splitting function of a parton j into a parton i which is 
related to the probability of a parton j emitting a parton % with longitudinal momentum 



3 



fraction z. These quark and gluon splitting functions in vacuum is evaluated by Altarelli 



and Parisi in 



111 ] at the leading order in coupling constant a s . 



In this paper we will evaluate the quark and gluon splitting functions in non-equilibrium 
QCD by using closed-time path integral formalism. We find that these splitting functions 
depend on non-equilibrium distribution functions of quarks and gluons in the QCD medium. 
The Q 2 evolution of the non-equilibrium fragmentation functions (eqs. ([1]) and ([21)) can be 
studied from eq. ([6]) by using these non-equilibrium splitting functions. 

We find the following expressions for the quark and gluon splitting functions in non- 
equilibrium QCD at leading order in coupling constant a s 

P gq {z) = C 2 {R) [1 + f q {k)f [1 + f g (k T , zk)f [1 + f q (-k T , (1 - z)k)f + Z)2 ] 

l + z\ 



p m (z) = (MR) [i + f q (k)Y [i + f g (-k T , (i - z)k)Y [i + f g (k T , zk)Y [- 

P gg (z) = 2C 2 (R) [1 + f g (k)} 2 [1 + f g (k T , zk)f [1 + f g (-k T , (1 - z)k)f 
1 ~* + T ^- + z(l-z)}. (7) 



z 

where 

C 2 {R) = \ (8) 

and k is the momentum of initial parton (which is assumed to be along longitudinal di- 
rection), kr is the transverse momentum of the emitted parton and z is the longitudinal 
momentum fraction of the initial parton carried by the emitted parton. 

Eq. ([7]) can be compared with the following expressions for the splitting functions in 



vacuum obtained by Altarelli and Parisi 

P gq {z) = C 2 (R) 



ll| at the leading order in coupling constant a & 
l + (l-z) 2 



P qq (z) = C 2 (R) ^ 

P gg (z) = 2C 2 (R) i 1 -^ l + T ^- z +^- *)]■ (9) 

Hence comparing with the splitting functions derived by Altarelli and Parisi in vacuum (see 
eq. (Q) we find that the non-equilibrium QCD splitting function (see eq. ([7j)) carries an 
extra factor [1 + fA(PA)] 2 x [1 + /b(pb)] 2 x [1 + fcipc)} 2 where parton A emits a parton 
B. 
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We will present derivation of eq. ([7]) in this paper. 

The paper is organized as follows. In section II we briefly review the derivation of quark 
and gluon splitting functions in vacuum. In section III we describe Schwinger-Keldysh 
closed-time path integral formalism in non-equilibrium QCD relevant to our calculation. In 
section IV we derive quark and gluon splitting functions in non-equilibrium QCD by 
using closed-time path integral formalism. Section V contains conclusions. 



II. QUARK AND GLUON SPLITTING FUNCTIONS IN VACUUM 

In this section we briefly review the derivation of quark and gluon splitting functions 
Pij in vacuum. We will present our calculation in the S— matrix approach. Hence, our 



derivation is slightly different from 111 ]. 



Consider a quark with momentum pa emitting a gluon with momentum ps in the process 
(lip a) —> g{pB) + live)- The S'-matrix element for this process is given by 



S« =ig J d 4 xN{i;ix)4 a ix)T a 4>ix)] (10) 
where (the normalization is from [13I ]) 



^(x) = ^ + (x) + $-{x) = EE \ T7^r\<P>ip)e- ip - x + at(p)u(p)e*- a 

spin P \ V -^P 



m 

spin /< v 



${z) = 4> + ix) + = J y^KipMp)^ + a\ip)uipy*- 



A»{x) = A^(x) + A»-{x) = EE J^rKW^)^ 1 + aj(p)^(p) e **]. (11) 

spin P I Z ^ fi P 

In the above equation a q (p), dq(p) and a g (p) are annihilation operators for quark, antiquark 
and gluons respectively. In eq. (fTTT) the suppression of color indices are understood. The 
initial and final states are 

\i >= \q{pa) >= a\( PA )\0 >, \f >= \q(pc), g(p B ) >= al(p o )al(p B )\0 >, (12) 

where Pc = Pa — Pb- 

By using eqs. GDI), and we find 



< f\S m \i >= '»/^E E^/^^ < OMpcKM 
[4(p)«(p)e i ''1[7„ot(p')£°''(p')e i '''-r-][a s ( P '')«( P ")e- i ''"-»] O ;^)|0 > . (13) 



For the x— integration in 

J d 4 xe ix ^ +p '- p "^ = J dt^ E+E '- E "^ J cPxe-^W-W, (14) 
we make a Wick rotation in time to find 



/ dW*W-"") = [ dte^ E+E '- E "^ ( d 3 xe~ i3 <^-^ = , - Z, f s (15) 
J Jt^oo Jv^oo (E + E'-E") v ; 

where we have used the fact that p+jf = p" in our partonic process. In the above assumption 
the error in our calculation goes to zero as the volume V — > oo 
The commutation relations are given by 

[a(p),aV)] = <W> [a( P ),a(p')} = [tf(p),tf(j/)) = 0. (16) 

Using eqs. (fT5|) and (TLB"]) we find from eq. ([13]) 



< f\S w \t>=tgJ2 E 



tr'„(E + E'- E") V U£„ V VE v n V 2V"E„/ 

spin p,p',p" \ ' / V P V P V P 

< 0|[<4(p)a 9 (p c ) + 8nv]u(p)jp[al(p')ag(p B ) + 5 plpB \ea^{p')T a [a\{p A )a q {p") + VpaMp")I° > • 



(17) 



Summing over p, p' and p" we obtain the transition matrix element square 

| < f\sW\i > | 2 = [- -l 2 — — — y |MI 2 . (18) 

where 

M = igu(pc)lMVA)e a ^p B )T a . (19) 

In order to obtain the probability for an initial quark emitting a final gluon in this process 
q(pa) 9{pb) + q(j>c), we need to multiply the phase-space factor of the remaining 

final state quark q(pc)- We find 

W 9q = | < 9(pbU(pc)\S<%(pa) > = I < f\S^\i > l 2 ^- (20) 



For quark spin sum we use 



E«(p)«(p) = ^-- (21) 

spin 
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Using eqs. QSJ and in p)} we find 



f22l 



For physical transverse gluons and massless quarks we find 

W~ - C ^^ 2E A E%E B?( E C + E B -E A? «« - >• (23) 

We assume that the quark p^ is along z— axis with the longitudinal momentum p for 
which we find 11] 



Pa — (p, 0, 0, p), p B = (zp + ^-,p T} zp) } pc = ((1 - z)p + ^ -p T , (1 - z)p) 

(pa-pbWc-pb) h , 2 P 2 T p _ . . 9,1 ^_ul£ ^ 
=2 = (1 - zjP - 7T^, ^A^C + Pa • Pc = 2(1 - z)p + — , (24) 

where z is the longitudinal momentum fraction of the initial parton carried by the emitted 
parton and pt is the transverse momentum of the emitted parton. Using this and performing 
the Dirac trace we find from eq. fl23l) 

u/ dzd Pr E B (1 - zf 2p 2 T 1 + (1 - zf f ^a s l + (l-z) 2 2 

W = ° 2{R)9 32^(1-^) 11 p\ z(l-z) z = ° 2{R) ^ z dZ d(lnPT) - 

(25) 



A. Quark to Gluon Splitting Function in Vacuum 

We find from eq. ( 1251) the quark to gluon splitting function 

Pgq ( z )=C 2 {R)—± L (26) 

which reproduces the exact result obtained by Altarelli and Parisi by a slightly different 
method [111 ]. 



B. Quark to Quark Splitting Function in Vacuum 

The quark to quark splitting function [in the process q{j>A) — ► q(pb) + gijPc)} can be 
obtained from the quark to gluon splitting function [in the process q(pA) 9{Pb) + q(pc)] 



7 



with the replacement z 




P qq (z)=P Gg (l-z), 



Z<1. 



(27) 



Hence we find from eq. (1261) the quark to quark splitting function in vacuum 



P„(z)=C*(R) 



1 - z 



(28) 



which reproduces the exact result obtained by Altarelli and Parisi. 
C. Gluon to Gluon Splitting Function in Vacuum 

Similarly, using three gluon diagram and carrying out the above algebra we find gluon to 
gluon splitting function in vacuum 



which reproduces the exact result obtained by Altarelli and Parisi by a slightly different 
method [111 ]. 

III. NON-EQUILIBRIUM QCD USING CLOSED-TIME PATH FORMALISM 

Unlike pp collisions, the ground state at RHIC and LHC heavy-ion collisions (due to the 
presence of a QCD medium at initial time t = ti n (say tj„=0) is not a vacuum state |0 > 
any more. We denote \in > as the initial state of the non-equilibrium QCD medium at t in . 
The non-equilibrium distribution function f(k) of a parton (quark or gluon), corresponding 
to such initial state is given by 




l-z 



z 



(29) 



< a ] (k)a(k') >=< in\a\k)a{k')\in >= 



(30) 



where we have assumed space translational invariance at initial time. 

Finite temperature field theory formulation is a special case of this when f(k) 




i 



S 



A. Quarks in non-Equilibrium 



The non-equilibrium (massless) quark propagator at initial time t = ti n is given by 
(suppression of color indices are understood) 

/ 



G(k)i 



k 2 + 



^ + 27r5(k 2 )f q (k) -27r5(k 2 )6(-k ) + 27r6(k 2 )f q (k) 



-2n5(k 2 )9(k Q ) + 27rS(k 2 )f q (k) 



k 2 - 



27rS(k 2 )f q (k) 



(31) 



where where i, j = +, — and f q {k) is the arbitrary non-equilibrium distribution function of 
quark. 



B. Gluons in Non-Equilibrium 



We work in the frozen ghost formalism 3j, |4j where the non-equilibrium gluon propagator 
at initial time t = t in is given by (the suppression of color indices are understood) 



[<T + (a - 1)— ] GJ^(k) - zT^G™ d (k) 



(32) 



where i, j = +, — . The transverse tensor is given by 

(k ■ u)(u^k v + u u k») - k»k v - k 2 u^u u 



V(k) = g 



(33) 



(k ■ u) 2 — k 2 

with the flow velocity of the medium . Gf^(k) are the usual vacuum propagators of the 
gluon 

(34) 



Gink) 



Wh, -2n5(k 2 )6(-k ) 

2\Q(1„ \ 1 



\-2nS{k 2 )B{ko) ^— 
and the medium part of the propagators are given by 



G™\k) = 2n5(k 2 )f g (k) 



1 1 
1 1 



(35) 



IV. QUARK AND GLUON SPLITTING FUNCTIONS IN NON-EQUILIBRIUM 
QCD 



In this section we evaluate quark and gluon splitting functions in non-equilibrium 



QCD. Similar to the vacuum case in 



ll| (see eq. ffl2]) ) we define the state \i > and |/ > in 
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non-equilibrium QCD as follows 

N >= \q(pa) >= a\{PA)\in >, \f >= \q(pc), 9(Pb) >= a\{pc)a) g {p B )\in > (36) 

where \in > is the initial state of the non-equilibrium QCD medium. It has to be remembered 
that for evaluating the Feynman diagrams and S— matrix we work in the interaction picture 
where the fields if>(x) and A^(x) obey the free field equations in terms of creation and 
annihilation operators as given by eq. (flip . From eqs. fl3"61) and ffTU]) we find 

< f\S w \t >= ig I d*x < tn\a q ( Pc )a g (p B )N$(x)4 a (x)T a ^(x)}al(p A )\tn > (37) 



which gives by using eq. ffTTT) 



< >= */AEEj^^ < M^c)a g [p B ) 

[at(p)«(p)e*l[7„4(pO£°"(pO^"T''][a„(p")«(p' , )e- i '''']at(p A )|m > . (38) 
Performing x— integration (see eq. ([To]) ) we find 



[aJ(p)«(p)][7X(? / )^(P') Tfl ][ a < r(P' , )«(p' , )]ai(PA)|tn > . (39) 



^„ (E + E' — E") V V£ p y V£> V 2 ^p' 
< :j \i>i" i i' i/'/-' j .<'./>/' i">// <•.'!'»» •- 

In the interaction picture the commutation relations are same as that for free field operators 

Hp), a\ P ')\ = <v, Hp), a( P ')\ = H( P ), a +(p')] = o. (40) 

which gives 



< f]sa)] * >= i9 ^(E + E--E») /^WITO < m| 

K(p) a «bc) + ^pcMPh^^'KM + Sp> PB ]e atl ( P ')T a [al( P A)a q {p") + 5 p>A ]u(p")\in > . 

(41) 

For our purpose of evaluating Feynman diagrams in momentum space we use eq. (130]) 

<in\ala p ,\in>=f(p) 8<§ (42) 

where we have assumed the space-translational invariance at initial time t — ti n — 0. Using 
eq. (142]) and summing over p, p' and p" we find from eq. (l4~Tj) 

1/ l9 m m 1 



</|S (1) |«> 



L (E c + E B - E A ) 1 VE C VE A 2VE B 
x [l + f q (pc)] 2 x [I + Upb]] 2 x [1 + /,(pa)] 2 £|M| 2 ( 4 3) 



spin 
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where 



M = igu(pchMPA)e a >*(PB)T a . (44) 

From now onwards we can follow exactly the same steps as in the vacuum case (see section 
II, the derivations after eq. (11811 ) to find the probability 

W = C 2 (R)^- [1 + f q (p c )} 2 x [l + f g (p B )} 2 x [l + f q (p A )} 2 \ ] dzd{W T ){Ah) 

Comparing this with the Altarelli-Parisi equation for the probability in vacuum (see eq. 
fl25|) ) we find that the probability in non-equilibrium QCD carries an extra factor [1 + 
Ia^Pa)] 2 x [1 + fsiPB)] 2 x [1 + fcipc)} 2 where parton A emits a parton B. 

A. Quark to Gluon Splitting Function in Non-Equilibrium QCD 

From eq. (j4"5j) we find the quark to gluon splitting function in non-equilibrium QCD at 
leading order in a s : 

P gq {z) = C 2 (R) [1 + f q {p)f [1 + f g (p T , zp)} 2 [1 + f q (-p T , (1 - z)p)f 1 + (1 ~^ )2 (46) 

which reproduces eq. ([7]). 

In the above equation p is the momentum of initial quark (which is assumed to be along 
longitudinal direction), px is the transverse momentum of the gluon and z is the longitudinal 
momentum fraction of the initial quark carried by the gluon. 

B. Quark to Quark Splitting Function in Non-Equilibrium QCD 

The quark to quark splitting function [in the process q{j>A) lips) + g(Pc)} can be 
obtained from the quark to gluon splitting function [in the process q(pA) 9(.Pb) + <?(Pc)] 
with the replacement z — > (1 — z) 

P qq (z)=P Gq (l-z), z<\. (47) 

Hence we find from eq. (1461) the quark to quark splitting function in non-equilibrium QCD 

P qg (z) = C 2 {R) [1 + f q {p)f [1 + f g (-p T , (1 - z)p)f [1 + f q (p T , zp)} 2 (48) 

which reproduces eq. ([7]). 

11 



C. Gluon to Gluon Splitting Function in Non-Equilibrium QCD 

Similarly, using three gluon vertex and carrying out the similar algebra we find gluon to 
gluon splitting function in non-equilibrium QCD 

P gg (z) = 2C 2 {R) [1 + f g {p)f [1 + f g (p T , zp)f [1 + f g (-p T , (1 - z)p)f [i^i + -JL- + Z (i- 

which reproduces eq. ([7]). 

Comparing eqs. (H6|) . (HHI) and (H9l) with the Altarelli-Parisi splitting functions in vacuum 
[see eqs. ( 1261) . f l28l) and ( 1291) ] we find that the non-equilibrium QCD splitting function carries 
an extra factor [1 + /a(pa)] 2 x [1 + }'b{.Pb)] 2 x [1 + /c*(pc)] 2 where parton A emits a 
parton 5. 



V. CONCLUSIONS 



RHIC and LHC heavy-ion colliders are the best facilities to study quark-gluon plasma 
in the laboratory Since two nuclei travel almost at speed of light, the QCD matter formed 
at RHIC and LHC may be in non-equilibrium. Since the fragmentation function is a non- 
perturbative quantity, we do not have theoretical tools in QCD to calculate it yet. The 
normal procedure at high energy pp, ep and e + e~ colliders is to extract it at some initial 
momentum scale /io and then evolve it to another scale \i by using the DGLAP evolution 
equation which involves splitting function Pji of a parton j into a parton i. These quark 



ll| at the 



and gluon splitting functions in vacuum is evaluated by Altarelli and Parisi in 
leading order in coupling constant a s . 

In this paper we have evaluated the quark and gluon splitting functions in non-equilibrium 
QCD by using closed-time path integral formalism. These splitting functions can be used 
to study DGLAP evolution equation in non-equilibrium QCD. For quarks and gluons with 
arbitrary non-equilibrium distribution functions f q (p) and f g (p), we have derived expressions 
for quark and gluon splitting functions in non-equilibrium QCD. We have found that the 
quark and gluon splitting functions depend on non-equilibrium distribution functions f q (p) 
and fg{p)- We have made a comparison of these splitting functions with that derived by 
Altarelli and Parisi in vacuum and have found that the non-equilibrium QCD splitting 
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function carries an extra factor [1 + /^(pa)] 2 X [1 + /b(pb)] 2 x [1 + fc(pc)} 2 where parton 
A emits a parton B. The Q 2 evolution of this non-equilibrium fragmentation function may 
be necessary to study hadron formation [8] from quark-gluon plasma Q,Q,Q,H at RHIC 
and LHC. 
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